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Abstract: The HARP system of resistive plate chambers (RPCs) was designed to perform par- 
ticle identification by the measurement of the difference in the time-of-flight of different particles. 
In previous papers an apparent discrepancy was shown between the response of the RPCs to min- 
imum ionizing pions and heavily ionizing protons. Using the kinematics of elastic scattering off 
a hydrogen target a controlled beam of low momentum recoil protons was directed onto the RPC 
chambers. With this method the trajectory and momentum, and hence the time-of-flight of the pro- 
tons can be precisely predicted without need for a measurement of momentum of the protons. It 
is demonstrated that the measurement of the time-of-arrival of particles by the thin gas-gap glass 
RPC system of the HARP experiment depends on the primary ionization deposited by the particle 
in the detector. 



- 1 - 



Keywords: Gaseous detectors, dE/dx detectors, particle identification methods, timing 
detectors . 



* Corresponding author, e-mail: Roumen.Tsenov@cem.ch. 



Contents 

|I[ Introduction ^ 

g HARP RPC system g 

|[ RPC time response g 

^ Discussion § 

|[ Conclusions ^ 



1. Introduction 

The HARP experiment has been designed to measure hadron production cross-sections on 

nuclear targets with a precision of a few percent over almost the full solid angle. A set of solid and 
liquid targets spanning a large range in atomic number was exposed to beams of protons and pions 
with momenta between 1.5 GeV/c and 15 GeV/c. The elements used ranged from hydrogen to 
lead. HARP took 450 milhon physics triggers, collected data for about 300 different settings and 
recorded 30 TB of information from August 2001 to October 2002. 

The setup of the HARP experiment is shown in figure [I| A detailed description of the detector 
and its performance is given in Ref. ||l|], [^. The spectrometer can be subdivided into three main 
systems. Beam and trigger detectors provide tracking and identification of beam particles, and 
trigger decisions. Forward detectors provide tracking, momentum measurement and identification 
of secondary particles at angles less than 17° with respect to the beam axis. Large angle detectors 
deal with tracking, momentum measurement and particle identification at large production angles. 

The large angle detector system consists of Time Projection Chamber (TPC), Resistive Plate 
Chambers (RPC) and an Inner Trigger Cylinder (ITC). 

The TPC has a cylindrical form with a length of ~ 200 cm and a diameter of ~ 83 cm. It is 
located inside a solenoid magnet producing a 0.7 T field. It measures momentum, trajectory and 
ionization energy losses of particles emitted from the target at large angles (20° < 6 < 160°) with 
respect to the incoming beam. The RPCs are arranged in the shape of a barrel around the outer 
field cage of the TPC. A target station and the ITC are situated inside TPC volume, in a truncated 
inner field cage (see figure ^. 

This paper concentrates on a specific effect found when calibrating the time response of the 
HARP RPCs. Besides the well-known time-walk correction, their time response was found to 
depend on the primary ionization deposited in the gas gap. Brief information needed to understand 
the RPC hardware and calibration procedures is given in the next section. The evidence for the 
effect itself is given in the Section |3[ followed by some discussion (Section 0) and conclusions 
(Section 1). 
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Drift Chambers 




Figure 1. Schematic layout of the HARP detector. The detector covers a total length of 13.5 m along the 
beam direction and has a maximum width of 6.5 m perpendicular to the beam. 




Figure 2. Schematic layout of the TPC. The beam enters from the left. Starting from the outside, first the 
return yoke of the magnet is seen, closed with an end-cap at the upstream end, and open at the downstream 
end. The cylindrical coils are drawn inside the yoke. The field cage is positioned inside this magnetic 
volume. The inner field cage is visible as a short cylinder entering from the left. The ITC trigger counter 
and the target holder are inserted in the iimer field cage. The RPCs (not drawn) are positioned between the 
outer field cage and the coil, as indicated. 

2. HARP RPC system 

Particle identification (PID) at large angles is mainly performed by mean energy loss {dE / dx) mea- 
surements in the TPC. A complementary system for particle identification to distinguish between 
electrons and pions in the momentum range 1 25 MeV / c-250 MeV /c is needed. For these relatively 
low momenta PID can be performed by time-of-flight (ToF) measurements. A time resolution of 
~ 200 ps is required when particles traverse the TPC over the shortest distance {r^ 400 mm). ToF 
measurements with such a resolution also help in pion-proton separation up to 1 GeV/c. A sys- 
tem of resistive plate chambers has been designed and used as additional PID device for particles 
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Figure 3. Cross-section of RPC glass stack. Six glass plates form four 0.3 mm gaps using fishing lines as 
spacer. The only signal electrode is situated in the middle of the stack. 



emerging with large production angles. 

The design of the HARP RPCs is based on a prototype developed for the ALICE experi- 
ment in, [§. 

The RPCs are constructed as a four-gap stack of glass plates. The gap size is precisely set 
to 0.3 mm by interposing a fishing line (nylon mono-filament) with suitable diameter between the 
plates. The stack consists of two identical structures of three glass plates each, arranged symmet- 
rically on both sides of a central readout electrode. The glass plates are 0.7 mm thick and made of 
standard float glass with a specific resistivity of ~ 10^^ Q.cm. A view of the cross-section through 
the short side of the glass stack is shown in figure ^ The negative high voltage is applied to the 
glass plates by means of a coated graphite layer with resistivity of 200 kD./\I\ on the two outer glass 
plates of both sets. A single readout electrode, located in the center of the glass stack, collects sig- 
nals from all four gaps. It is segmented into 64 rectangular strips. Eight strips are connected to 
one pre-amplifier forming a readout channel. Thus, each chamber has eight readout channels, in 
the following referred to as pads, and numbered from 1 to 8 starting from the most upstream one. 
Pads with one and the same number define a pad-ring. A chamber is 1920 mm long, 104 mm wide 
and 7.8 mm thick. Each chamber is housed in an aluminum box with dimensions 2 m x 150 mm 
X 10 mm. 

Thirty such chambers are arranged around the TPC in two staggered layers forming a barrel. 
They cover polar angles from 17° to 142° with respect to the beam axis and 2n in azimuthal angle 
with a small 13 mm overlap between the layers. The set of two layers fits into a ~ 25 mm radial 
space between the TPC and the coils of the solenoid magnet. The readout electrodes of the inner 
and outer layer are located at radial distances of 421 mm and 434.5 mm, respectively. The front- 
end electronics of each channel consists of a pre-amplifier board mounted on the chamber and a 
combined splitter/discriminator module. 

Pre-amplifiers are based on the AD8009 chip operating with an amplification factor of ~ 30. 
The amplified signals are transmitted through mini coaxial cables over a distance of 0.8 m - 2.5 m 
(depending on the channel in question) to a passive patch panel and from there over a distance of 
5 m through Lemo 50 D. cables to a custom-made splitter and leading edge discriminator module. 
Each signal is amplified once more there (with a factor ~ 3 using the same AD8009 chip) and split 
into two separate signals. One of them is discriminated (the discriminator threshold is put at 5 mV) 
and is sent via a 80 m twisted pair cable to a Time-to-Digital Converter (TDC), model CAEN 
V775 with a channel nominal width of 35 ps. The other signal is sent through another 80 m twisted 
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pair cable to a Charge-To-Digital Converter (QDC), model CAEN V792. 

The RPCs were operated in avalanche mode at a voltage of —6 kV between outer and central 
electrodes and with a gas mixture of 90% C2F4H2, 5% SFg, and 5% C4H10. 

Typical random noise rates were in the range of 200 Hz-300 Hz per chamber, i.e. < 1 kHz/m^, 
which is an acceptable level compared to the typical particle rate 10 kHz/m^ over the area covered 
by the barrel RPC. 

A precise calibration of the RPC sub-detector is needed before it can be used for time-of-flight 
measurements. The aim of the calibration is to develop a procedure that transforms data read out 
from a TDC channel into a flight-time of particles from the production point in the target to the RPC 
pad feeding that particular TDC channel. An in-situ calibration procedure has been developed [||] 
by using reconstructed charged particle tracks in the TPC from interactions produced by beams of 
positive particles with momenta of 3 GeV/c, 5 GeV/c, and 8 GeV/c bombarding 0.05 Ajjjt (nuclear 
interaction length) thick Ta, Pb, Sn and Cu targets. Pion tracks were selected by using dEjdx 
measurement in the TPC. The algorithm converts the TDC scale into physical units, accounts for 
the arrival time of the beam particle at the target, for the transit time of the signal within the pad 
and from the pad to the preamplifier, and for temperature fluctuations of the time response. It also 
corrects for so called time-walk arising from the fact that signals with different pulse-height and 
the same shape cross the fixed discriminator level at different time. It is worth mentioning that 
the calibration procedure corrects for threshold crossing delays under the assumption of a single, 
universal pulse-shape of the signals from all particles. The correction is quite large, up to 2 ns and 
is of the order of the smallest time-of-flight to be measured. 

More details of the design, layout, calibration and operational parameters of the HARP RPCs 
can be found in [f], [0], %. 

3. RPC time response 

Figure ^ demonstrates the correlation between measured relativistic velocity jS of positively charged 
particles and their momentum measured in the TPC when the time-walk is calibrated with pion sig- 
nals and applied to all tracks It is clearly seen that proton signals are shifted with respect to the 
theoretical curve. The rest of this section is devoted to the understanding of the origin of this shift. 

In Ref. [^] two mechanisms were suggested which could be the cause of this effect. One 
possible explanation is the fluctuation in arrival time of the first cluster of the primary ionization. 
This fluctuation is smaller for heavily ionizing particles. The other possibility is a change of pulse 
shape near threshold due to a possibly different ionization and different gas amplification regime. 

Since the observed rise-time of the signal to the discriminator threshold is several ns, a small 
difference in signal shape between heavily ionizing particles (protons with momenta between 
250 MeV/c and 800 MeV/c in our case) and minimum ionizing particles (here pions with mo- 
menta above 150 MeV/c) could induce such an effect. Effects below a ns are contained within the 
expected avalanche formation time. 

Dedicated analyses have been performed to understand the effect better. 

First, we compared the relative time offset between the measured over-threshold time of the 
signal and the predicted arrival time of protons above 1 GeV/c with relativistic pions and with 
heavily ionizing protons as a function of the measured energy loss of particles in the TPC gas 
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Figure 4. Relativistic velocity /3, measured by the barrel ToF system, as a function of momentum, measured 
by the TPC, for positive particles. The curves represent theoretical dependence for pions (upper curve), 
kaons (middle curve) and protons (bottom curve). Proton signals are shifted towards higher j3. 



(Figure To reduce the possible effect of the time-walk correction we use a total pulse charge 
range 1400 < Q < 1600 in terms of channel counts, which is a region of relatively high pulse 
charge where still a large statistics can be obtained. The energy loss in the TPC gas is correlated 
with the ionization density of the same particle crossing the RPC gas gaps, just shifted by the 
material between the TPC gas volume and the RPC gas gaps. Figure || shows a clear dependence 
of the time offset on the ionization density. The points where both pions and protons are minimum 
ionizing show a similar time offset, while the points at higher dE/dx display a trend justifying our 
hypothesis. 

This result came as first evidence that a dependence on the ionization density in the gas gaps 
plays a role |^]. 

Further, if our hypothesis is correct, the time difference between the observed signal (corrected 
for the time-walk per pad) and the expected arrival time calculated on the basis of the track length 
and particle momentum as measured in the TPC should show a dependence as a function of mo- 
mentum for the same set of pad-rings due to the correlation between the energy loss and momentum 
for a given particle species ' . Such a dependence was suggested in Ref. [||] and the respective plots 
are reproduced here in figure ^. The behavior seen in the figure indeed demonstrates a clear de- 
pendence on ionization density (dE/dx) and a dependence on different path lengths at fixed dE/dx 
only for low values of this quantity. 

It should be noted that the above momentum dependencies were obtained using the momentum 
of the particles measured in the TPC gas. In Ref. [[I0|], | U] it has been pointed out that a similar 



behavior can be obtained when a systematic shift in the measurement of momentum is present. 



'it is important to distinguish between different pad-rings because particles emerging from the target cross the pad- 
rings at different inclinations, hence they produce different total ionization even having the same momentum. 
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Figure 5. Offset of the time-over-threshold of RPC signals selected with measured total charge between 
1400 and 1600 channel counts as a function of dE /dx in the TPC gas. The dE jdx scale is in arbitrary units 
with the value of the minimum ionizing particle at « 150 counts. The closed circles are obtained with pions, 
the open circles with protons using ?l dE /dx selection, and the open square using protons with momentum 
above 1 .5 GeV /c in addition to the dE / dx selection. The points where both pions and protons are minimum 
ionizing show a similar time offset. 
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Figure 6. Difference of the measured time offset from the expected time offset for pions (left panel) and 
protons (right panel) as function of the momentum measured in the TPC, for pad-rings 3, 5, and 7. 



Besides our confidence, based on various independent calibration methods, that our momentum 
measurement is unbiased (||^, [^, [|l2|], [[l^]) the effect of any momentum measurement bias in 
the TPC can be eliminated by an analysis employing the kinematics of elastic scattering using a 



-6- 



liquid hydrogen target. Such a measurement makes it possible to send a "controlled beam" of slow 
protons through the TPC and towards the RPC system without the need to measure the momentum 
of the recoil proton with the TPC. The prediction for the momentum and direction of the recoil 
proton can be obtained from the kinematics of the event by measuring the scattering angle of the 
forward scattered proton or pion. Exploiting this, we have used an exposure of the HARP detector 
where a 5 GeV/c beam of protons and pions is directed onto a 60 mm long liquid hydrogen target. 
Elastic events are selected from the total sample of triggers requiring one and only one track in 
the forward spectrometer, and exactly one track in the TPC. The momentum of the track measured 
in the forward dipole spectrometer is required to be consistent with elastic scattering hypothesis. 
Only events with exactly one RPC hit are retained. The position of the RPC hit has to be consistent 
with the predicted impact point of the recoil track, using the direction of the forward scattered 
track to define the trajectory of the recoil particle. This selection of elastic scattering events has a 



purity of 99% [|12|]. Due to the acceptance of the forward spectrometer, the distribution of selected 
recoil protons peaks at 74° with respect to the beam direction, thus almost perpendicular to the 
RPC chambers, with most of the tracks contained between 60° and 80°. The recoil momenta are 
between 350 MeV/c and 600 MeV/c. The geometry of the RPC system is such that the large 
majority of selected recoil protons traverses pad-ring 3, therefore, this measurement is performed 
using this pad-ring only. In order to predict the time-of-flight of the particle reliably, the initial 
momentum and energy loss in the material traversed has to be described accurately. For example, 
a proton with 325 MeV/c momentum in the TPC gas has lost in terms of momentum on average 
~ 45 MeV/c with an r.m.s. of 15 MeV/c in the material between the interaction point and the TPC 
gas volume and is going to lose 83 MeV/c more (with an r.m.s. of 25 MeV/c) in the outer field 
cage before reaching the RPC detection volume. At 575 MeV / c the total reduction in momentum 
is ~ 38 MeV/c on average with similar r.m.s. ^ 

The properties of the kinematics of elastic scattering are used to predict the path-length of the 
particle from the interaction point to the RPC detector and the time-of-flight over this distance. 
Monte Carlo simulation is used to verify that the predicted and reconstructed time and path-length 
agree within 5 ps and 2 mm, respectively. In the left panel of figure ^ the time difference, defined as 
the measured time minus the predicted time, is displayed as a function of the predicted momentum 
of the proton in the gas of the TPC. The momentum was predicted using the kinematics of elastic 
scattering. The fact that the simulated difference of prediction and measurement is consistent with 
zero shows that the prediction of the flight-time (and thus of the momenta) using the elastic scatter- 
ing kinematics and Monte Carlo corrections in the reconstruction procedure for respective energy 
losses is correct. The data exhibits a clear deviation pointing to a difference in RPC time response 
as a function of the momentum. The right panel of figure shows the same time difference as a 
function of the momentum predicted at the RPC detection volume. Plots in figure ^ clearly demon- 
strate the existence of a remarkable dependence of the RPCs response on the primary ionization 
density produced in the gas-gap by the particles of different types. The magnitude of the effect is 
in agreement with our previous analysis [^, [||], [^]. 

In fact, one would expect that the RPC response does not depend directly on momenta of 



^The general HARP Monte Carlo simulation package (see ||2||) based on GEANT4 toolkit [0] has been used through- 
out the analysis presented in this article. 
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Figure 7. Left panel: The difference of the time offset measured in pad-ring 3 from the expected time offset 
for protons as a function of the momentum along its flight path (in the gas volume of the TPC). The filled 
circles show the results of measurements using elastic scattering on hydrogen, the points without marker 
represent the simulation of the measurement using the same reconstruction procedure. The momentum was 
predicted using the kinematics of elastic scattering. Consistency of the simulated time difference with zero 
shows that the prediction of the flight time (and thus of the momenta) using the elastic scattering kinematics 
and Monte Carlo corrections in the reconstruction procedure for respective energy losses are correct. 
Right panel: The difference of the time offset measured in pad-ring 3 from the expected time offset for 
protons as a function of the momentum in the RPC gap. 



detected particles, but only on their energy deposition in the gap. Given the small range of path 
lengths in the RPC gas gaps in the sample of elastic events it is natural to show the time differences 
as a function of total energy deposition in the RPC gas. The energy deposition can be predicted as 
a function of momentum using our simulation. In figure ^ the same time difference measurements 
as the ones of figure ^ are displayed as a function of the above energy deposition. The range of this 
quantity in the data is between two and eight MIP (the average deposition of one minimum ionizing 
particle). In this limited range a nearly linear behavior of the shift in time response is observed, with 
the extrapolation to one MIP consistent with a vanishing shift. This is to be expected, since the RPC 
system was calibrated using pions with momenta close to the value where their energy deposition 
is minimal. Although the dependence does not show a flattening at higher energy deposition in the 
measured range, one would expect this to occur at very high ionization losses. 



4. Discussion 



We have clearly shown in the previous section that the time response of RPCs with thin gas gaps 
depends on ionization density of detected particles in those gaps even after correction for the usual 
time-walk based on the total charge of the signal. Evidence for such an effect was reported by 
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Figure 8. The difference of the time offset measured in pad-ring 3 from the expected time offset for elasti- 
cally produced recoil protons as a function of the predicted energy deposition in the RPC gaps. The energy 
deposition is expressed in units of the average energy deposition of a minimum ionizing particle (MIP). 



us earlier [g], [15]. Now we support our findings by strong additional evidence exploiting 

proton-proton and pion-proton elastic scattering kinematics that allows us to avoid any dependence 
on the momentum measurement in the TPC. 

In view of the evidence presented in this work, the comments of Ammosov et al. about our 



first indication M\ of the effect, e.g. [10|: 



The timeslewing correction for the timing of pulses with finite rise time has been 
well understood for over half a century. The discovery that on top of the timeslew- 
ing correction a further correction for protons 'due to their higher ionization rate' is 
needed, is only one example out of many grave misconceptions and mistakes contained 
in the above-cited paper. ... The above-mentioned discovery that the proton timeslew- 
ing correction is different from the pion timeslewing correction by 500 ps is a plain 
mistake stemming from the lack of understanding of the interplay between TPC track 
and RPC timing reconstruction. 

show that for the authors of [ |To| ] such an effect was unexpected and therefore missed. 

It is interesting to mention that, after more evidence for the discussed effect was published in 
our papers [Ql, [^, the same authors acknowledged the effect without proper reference to our work 



and tried to construct a physical model explaining it [ 16]. Their estimation of 80 ps time difference 



for protons at 400 MeV/c is by factor of 2.5 smaller than our measurement. In a very recent 



preprint [ ]17p the same authors exploit the time-of-fiight of low momentum protons measured by 
the same RPC system as a benchmark for their procedure of momentum measurement in the HARP 
TPC. They do not account for the effect reported here risking a shift in their momentum calibration. 
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5. Conclusions 



It was found that the time response of the thin gap RPC system used in the HARP experiment 
and caUbrated with the assumption of a unique pulse-shape appeared to be different for pions and 
protons in the momentum range below 1 GeV/c. In earlier publications we proposed the higher 
energy loss of protons as a probable explanation of the effect. In this paper it was clearly confirmed 
that the effect is present and can indeed be understood as originating from the considerably higher 
ionization produced by the protons in the RPC. Momentum measurement biases in the TPC, if any, 
have been ehminated as possible cause of the effect. 

It should be noted that possible measurement of the effect using a direct beam of protons with 
momenta as low as 300 MeV/c-600 MeV/c would be interesting. However, it is not an easy task 
due to the high energy losses of such protons. 
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